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We experimentally demonstrate the broad tunability of the main features of optical 
' localized structures (LS) in a nonlinear interferometer. By discussing how a single 
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LS depends on the system spatial frequency bandwidth, we show that a modification 
of its tail leads to the possibility of tuning the interactions between LS pairs, and 



CN ' thus the equilibrium distances at which LS bound states form. This is in agreement 

> : 

with a general theoretical model describing weak interactions of LS in nonlinear 



dissipative systems. 

Localization of spatial patterns is a subject of major current interest in the re- 
search on nonlinear dissipative dynamical systems. The studies about this topics have 
naturally followed and sided those dedicated to the formation of temporal and spatial 
solitons in Hamiltonian systems Analytical and numerical works have identified 
several distinct mechanisms leading to structure localization in dissipative systems 
0, and experimental observations of this phenomenon have been recently offered in 
several systems, such as fluid dynamics [Q], chemistry Q, granular materials || and 
nonlinear optics 

In particular, optical localized structures (LS), to which we will also refer to as 
dissipative solitons in the following, are objects of intense research, also in view of 
possible applications as pixels in devices for information storage or processing. So 
far, the existence of optical dissipative solitons has been theoretically predicted in 
many passive @ and active || configurations, and optical LS have been observed in 
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photorefractive cavities || and in passive nonlinear interferometers, based either on 
the "thin slice with feedback" scheme [H], |11], |l2fl , or on a microresonator filled with 



a semiconductor medium [13|. More recently, the interactions between LS have been 



shown to give rise to the formation of a discrete set of bound states |TT 



To our knowledge, very little is known about the dependence of the LS's features on 
the experimental parameters. The present work addresses this issue, by investigating 
how the spatial frequency bandwidth of a nonlinear interferometer can be utilized to 
tune both the spatial profile of each single soliton, and the interaction forces occurring 
between two of them. A quantitative experimental evidence is given of the crucial 
role played by the oscillatory tails of a single LS in determining the interaction forces 
between solitons. 

Our experimental system consists of a Liquid Crystal Light Valve (LCLV) closed 
in an optical feedback containing both interferential and diffractive processes. When 
an initially plane wave is sent into the system, its phase (f(r,t) evolves according to 
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^ = ~(<P ~ Vo) + Ifrlv + otl | e^T (Be- 1 * + C) | 2 (1) 
where tpo = it is the phase working point of the LCLV, and r and Id are its response 
time and diffusion length respectively. The source term in the right hand side of Eq. 
(HD depends on the free propagation length I in the feedback loop, as well as on the 
laser light wavenumber k and on the parameters B and C, that tune the relative 
weight of diffraction and interference in the system. Finally, Iq is the incident laser 
intensity, and a describes the Kerr-like response of the LCLV. Here, B = cos 2 9, 
C = sin 2 9, where 9 is the (experimentally adjustable) angle between the director 
of the nematic liquid crystals of the LCLV and the transmissive axis of a polarizer 
oriented along the polarization direction of the incident light. 



In a previous work ||12|| , we have characterized the state diagram of the interfer- 
ometer in the parameter plane (9, Iq), finding that localization of patterns occurs for 
a broad range of 9 values (~ 35° to 58°). This phenomenon is related to the presence 
of a subcritical bifurcation, connecting a lower uniform branch to an upper patterned 
one. In these conditions, the formation of isolated spots connecting the two branches 



is typical |7|, [14], [15| . Besides 9 and I , the scenario of observable patterns crucially 
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depends on the spatial frequency bandwidth q B of the interferometer, which can be 
experimentally controlled by means of a variable aperture put in a Fourier plane. In 
what follows we discuss the main LS features that emerge by keeping fixed 6 = 42°, 
and varying Iq and the adimensional parameter q b = qs/qdm obtained by normalizing 
the system bandwidth to the diffract ive interferometer wavenumber gdiff = \J ^f^- 

A first point of interest is to establish the range of existence of LS in the (<#,, Iq) 
plane. In Fig. 1 we plot the state diagram of the system in this parameter plane, 
together with some snapshots representative of the observed patterns. All the exper- 
iments are performed at incident laser wavelength A = 632 nm and for / = 250 mm. 
This results in a scale of the observed patterns of the order of 27r/gdifr — 0.5 mm. 

Looking at Fig.l, one easily realizes that the range of existence of LS is very broad, 
not at all limited to some particular parameter choices. The lower threshold for the 
existence of LS increases for decreasing q b . This is a consequence of the fact that LS 
have an internal structure conatining both low and high frequency components, as it 
will appear evident in the following. Therefore, any bandwidth limitation perturbs 
the LS structure, and increases the threshold for their existence. At very low qb and 
high intensities, localization of structures is lost and regular hexagons are observed, 
due to the long range correlation imposed to the pattern by the small bandwidth. 

If Iq is kept fixed at high values while q b is increased, hexagonal patterns evolve into 
a space-time chaotic (STC) regime. The boundary line between STC and LS occurs 
at decreasing intensities when q b is increased. This indicates that the regime here 
generically referred to as STC can arise either from a strong excitation of a relatively 
small band of wavenumbers, or from a weak excitation of a large set of interacting 
spatial modes. The indetermination of the boundaries between the different regimes 
is of the order of 10 %. It must be also specified that the placement of the boundaries 
depends on the evolutionary history of the parameters, since we are in presence of a 
subcritical bifurcation. The continuous lines in Fig. 1 were obtained by decreasing 
the input intensity, the dashed line by increasing it. Localized structures are not 
observed in this last case. 

Scanning the parameters within the domain of LS's existence leads to sensible 
modifications in the shape of each structure. In Fig. 2 we show the variation in 
the LS intensity profile observed by keeping Iq close to the lower threshold for LS 
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existence and increasing q b . It is seen here that each structure is formed by a central 
peak, and by a set of concentric rings forming a tail that shows spatial oscillations of 
decreasing amplitudes for increasing distances from the LS center. The width of the 
central peak can be roughly evaluated as the diameter of the first dark ring in each 
frame, and appears to be practically independent on q b . 

The length scale of the oscillations on the tails is instead strongly dependent on 
q b . Namely, this scale decreases for increasing q b until q b — 3, and then saturates to 
a constant value. 

The set of our observations indicates that LS have a "natural" unperturbed shape 
like that displayed for q b > 3. By constraining the system to a bandwidth smaller 
than this value, one is then able to tune the LS profile, imposing oscillations on the 
tails at a frequency different from the natural one. The occurrence of oscillatory tails 



on LS have been reported in other physical systems fl5| , \L% , and it is considered to 
be a typical signature of the formation of LS via pinning of the fronts connecting the 
uniform and the patterned states @. 

The observed LS closely resemble those reported in Ref. [HJ, in which a subcritical 
real Swift-Hohenberg (S-H) equation is studied analytically and numerically. This is 
not surprising, since our experiment displays a subcritical bifurcation of a real order 
parameter to a patterned state, and therefore is appropriately modeled by an order 
parameter equation of that kind. We do not expect that the S-H model describes 



faithfully all the details observed in the experiiment, however. It is known fig] , for 
example, that the "thin slice with feedback" model, of which our experiment is an im- 
plementation, presents instabilities at multiple wavenumbers given by = V^q&fi, 
N = 1,5,9,.... Though the highest wavenumbers become active at high values of 
pump parameter due to diffusion, it may be expected that they play some role in 
determining the fine features of the LS's. Our aim in comparing the experimental 
findings wit the prediction of the S-H model is indeed to investigate wheter some 
fundamental features of the observed phenomena can be described in terms of this 
very general model. 

Using the Swift-Hohenberg model, it is found analytically that the LS tails are 
described by single spatial scale oscillations, embedded in an exponential envelope 
that departs from the lower uniform state. 
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Though the LS tails in our case display some deviations from the above ideal 
behavior, the qualitative agreement between our observations and the results of the 



general theory reported in Ref. |T5[ is satisfactory. In particular, it is possible to 
identify for each value of a dominating spatial scale in the oscillatory tails. To 
this purpose, we measure the distance between successive maxima of a single LS and 
average this quantity over all observed maxima. This way, we obtain the dominant 
spatial frequency of the tail oscillations, which is then normalized to gdifr and reported 
as g ta ii s in Fig. 3. The error bars correspond to the measured frequency fluctuations 
from the gtails, reflecting the fact that the tail oscillations are not rigorously at a single 
spatial scale. Looking at Fig. 3, one easily realizes that g ta iis practically coincides 
with qb for qb < 3. At higher values of qb, no variations in q ta n s as well as in overall 
LS's profile are observed. 

The shape of the tails is responsible for the interactions between localized struc- 
tures. Namely, while for monotonically decreasing tails, one would expect only at- 
tractive or repulsive forces between LS, oscillatory tails induce oscillatory signs of the 
interactions, thus producing both attractive and repulsive forces, depending on the 
distance between the centers of a pair of LS's [|1|, [TTj] . A recent work [|11|] has exper- 
imentally demonstrated the existence of a discrete set of LS bound states, occurring 
in the presence of oscillations on the LS tails. The selection rule for the discrete set of 
bound states observed has been there put in relation with the spacing of their rings 
originated by diffraction around the central peak. In the following we show how these 
bound states can be in fact tuned by varying the spatial frequency bandwidth of the 
interferometer, and we discuss how the selection of the observed bound ststes can be 
put in the very general framework of a subcritical S-H model. 

In Fig. 4 we display a set of different bound states observed for qb = 3.6. We no- 
tice that the states form a set that can be ordered following a precise rule, given by 
simply counting of the number of maxima and minima that occur along the segment 
connecting the two LS centers. We will call this number n as bound state order num- 
ber. Such a feature is encountered for all values of qb- At small system bandwidths, 
however, we observe only the first two or three bound states, instead of the entire set 
shown in Fig. 4. This is probably due to the fact that the binding energy of each 
state varies with qb, and in some cases it is not sufficient to keep the LS pair tightly 
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bound in the presence of unavoidable system inhomogeneities and fluctuations. 



A theory for the interaction of LS pairs was given in Ref. [15[] in the context 
of study on a real Swift-Hohenberg equation. As already discussed, we expect this 
model to be closely applicable to our system in the present conditions. Following 
that approach, the weak interactions between a pair of LS with tails decaying with a 
length \i and oscillating at a frequency z/, lead to a time evolution for the distance R 
between the LS centers ruled by 

^ = ^ (e "" Rc ° s(z/jR)) = F{ ~ R ^^- ( 2 ) 
As a consequence, an infinite number of stable bound states are possible, corre- 
sponding to the solutions F — 0, dF/dR < of Eq. (2). In the limit in which the 
scales n and v are well separated, the difference R n +i — R n between the separation 
distances of two successive bound states corresponds approximately to the tail oscil- 
lation length v~ x of a single LS. We recall that by weak interaction we mean a regime 
in which the intensity amplitude of one LS is small in the space region in which the 
intensity amplitude of the other is large. In the case of our experimental data, this 
is true for all the bound states observed, with the possible exception of the lowest 
order one. 

If we assume that Eq. (2) describes correctly the bound state selection rule in our 
experiment, it immediately follows that tuning of the equilibrium distances should 
be possible by varying the scale of the oscillations on the tails of each single LS. 
In order to check this point, we measured the quantities A„ jn+ i = R n +i — R n , and 
then averaged them over the bound state order number n. The resulting quantity 
d (normalized to the length A ta ii s ) is reported vs. qb in Fig. 5. A constant value of 
the ratio -r-^ - — 1 is observed within the errors, indicating that the above discussed 

Ataiis ° 

relation between the oscillations on the tails of each LS and the selection rule of 
bound states is verified. This marks the fact that tuning of the equilibrium distances 
between LS's in bound states can be quantitatively performed in our experiment. 

In conclusion, we have given a quantitative evidence of the tuning of the LS spatial 
profile in a nonlinear optical interferometer, using the system spatial frequency band- 
width as a control parameter. We have discussed the role of the oscillations occurring 
on each single LS tail in determining the interactions between different LS's. Finally, 
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we have verified the agreement between the selection rules for the formation of bound 
states observed in our experiment, and those predicted for the same phenomenon by 
a general model for pattern formation in nonequilibrium systems. 
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FIG. 1: State diagram of the system in the ((#,, Iq) parameter plane. Notice that LS emerge 
for a broad range of parameters (gray area in the plane). The three reported patterns are 
snapshots of the observed hexagons, LS, and space time chaotic states. 
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FIG. 2: Variation of the LS's shape with the system bandwidth. Snapshots of the observed 
solitons for (a) q b = 1.0, J = 700 fiW/cm 2 (b) q b = 1.2, J = 620 fiW/cm 2 , (c) % = 
1.6, 7 = 520 fiW/cm 2 , (d) q b = 2.2, J = 480 ^VF/cm 2 , (e) % = 2.8, J = 460 ^VF/cm 2 , (f) 
q b = 4.0, 1 = 460 /iW/cm 2 . 
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FIG. 3: Variation in the main frequency of the LS tail oscillations as a function of the 
system bandwidth. Both g ta ii s and qb are adimensional quantities (see text for definition). 
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FIG. 4: Snapshots of different bound states observed at qj, = 3.6, Jo = 500 ixW/cm 2 . All 
patterns (a-e) are obtained by inducing a pair of LS with an increasing initial distance 
between centers, and letting the system evolve up to the time at which the stationary 
bound state is realized. 
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FIG. 5: 



vs. qb (see text for definitions). Both quantities are adimensional. Notice 



that, for all measurement a constant value of At - ~ 1 is realized within the experimental 
errors. 



